Half-vortices in polariton condensates 



Yuri G. RubcE 

Centro de Investigation en Energia, Universidad National Autonoma de Mexico, Temixco, Morelos 62580, Mexico 

(Dated: May 4, 2007) 

It is shown that vortices in linearly polarized polariton condensates in planar semiconductor 
microcavities carry two winding numbers (k, m). These numbers can be either integer or half-integer 
simultaneously. Four half-integer vortices (1/2, 1/2), (-1/2, -1/2), (1/2, -1/2), and (-1/2, 1/2) are 
anisotropic, possess the smallest energy, and define the Kosterlitz-Thouless transition temperature. 
The condensate concentration remains finite within the core of half-vortex and the condensate 
polarization becomes fully circular in the core center. 

PACS numbers: 71.36.+C, 42.55.Sa, 03.75.Mn 



Introduction. — Recent observations of exciton- 
polariton condensation in semiconductor microcavities 
U [|J revealed the formation of condensates with a 
well-defined linear polarization. The polarization is 
build-up as a result of minimization of the energy of 
polariton-polariton repulsion H lnt 0, 0| , 

Hint = I y'd 2 r{(C/o-C/i)(t/>*-^) 2 + {/ 1 |t/>*x^| 2 }. (1) 

Here the integration is over the microcavity plane within 
the excitation spot and the polariton condensate wave 
function (the order parameter) is written as a complex 
two-dimensional vector ip(r). This vector describes the 
in-plane component of electric field of the polariton con- 
densate and it is normalized to the condensate concen- 
tration n — (xp* ■ tp). The polariton repulsion is^ charac- 
terized by two interaction constants, Uq and U\ 4]. They 
are typically related as Uq > U\ > 0, so that at a fixed 
polariton concentration the minimum of iJi nt is reached 
at tp* x ip — 0, i.e., at a linear polarization. 

While the condensation has been observed in P, 0| 
for the case of localized polaritons, it is of key impor- 
tance to understand the structure and polarization prop- 
erties of topological defects (vortices) in uniform polari- 
ton condensates. First, it is vortex-antivortex unbinding 
that defines the critical temperature T c for the condensa- 
tion transition in two-dimensions (2D) according to the 
Kosterlitz-Thouless scenario [f|. And the estimations of 
T c given till now [|| did not take into account the polar- 
ization degree of freedom of polaritons. Secondly, since 
vortices are topologically stable objects, they can be used 
as long-living optical memory elements. The vortices 
studied below can be applied for polarization sensitive 
optical computing. 

Winding numbers and vortex interactions. — I begin 
with the classification of vortices in polariton conden- 
sates. Following the general approach [7j, it is necessary 
to distinguish topologically different mappings of a closed 
path in the order parameter space into a closed path 
around the vortex core in the microcavity plane. Since a 
linear polarization vector can be written in general as 

■0iin = {^x,i> y } = %/ne' e {cos?7,sinf7}, (2) 



the order parameter is defined by two angles, r](r) and 
0(r). The topology of the order parameter space, how- 
ever, does not coincide with the topology of a torus, be- 
cause the points 77, 9 and 77 + tt, 9 + n should be identified. 
It is more convenient to imagine the order parameter 
space as an infinite chessboard shown in Fig. [1] A vortex 
is defined by a line connecting two equivalent points on 
this chessboard, and it is characterized by two winding 
numbers, k and m, such that the angles are changed as 
t] — > T] + 2nk, 9^9 + 2irm by encircling the vortex core 
around. It is seen that the winding numbers k and m can 
take either integer or half-integer values together. Basic 
vortices carry half- integer values k,m = ±1/2. Corre- 
sponding changes of the order parameter are shown by 
four thick arrows in Fig. [1] Other vortices are reduced 
to superpositions of four basic half- vortices. 
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FIG. 1: Illustrating the order parameter space of polariton 
condensate. Blue arrows show Re{i/>}. Similar points in the 
black or white squares (as, e.g., the points marked with small 
crosses) represent identical values of the order parameter. 
Thick red arrows show four topologically different changes 
of the order parameter when one encircles the cores of four 
basic half- vortices. 
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Existence of half-integer vortices is a general feature of 
multi-component condensates @ , but usually they carry 
only one half-integer topological charge originating from 
the phase angle change. E.g., the half- vortices in 3 He- 
A generate rotations d — ► — d of the spin quantization 
axis d, accompanied by the phase change 9 — > 9 ± it @. 
Similar half-vortices have been discussed for the three- 
component spinor s = 1 atomic condensates [loj ]. For 
a 3D vector d, different trajectories connecting two an- 
tipodes on the sphere, d and — d, can be continuously 
transformed one to another, so that all half-vortices with 
the phase changing by 7r are topologically equivalent. 
Contrary, the polariton condensates are two-component, 
and the clockwise and counterclockwise rotations of the 
2D real vector {cos 77, sin 77} are topologically distinct. 
This gives rise to the second topological charge. 

Half-vortices with two winding numbers have b eing 



predicted for 3 He-A in the parallel-plate geometry [11 1 
where the spin quantization axis d is restricted to be 
two-dimensional. In this case, however, the sign of phase 
change is fixed and only two half-vortices can be present 
in the same configuration, while all four half-vortices can 
coexist in the polariton condensate. 

The total energy of polariton condensate is 



H 



2m 



Hint, (3) 



where m* is the polariton effective mass at the bottom 
of the lower polariton branch and p = (Uq — Ui)n [4[ 
is the chemical potential. The vortex energy consists of 
the core energy and the elastic energy E c \. The latter 
comes from polariton kinetic energy far away from the 
core, where Eq. © can be used. This gives 



~ 2 Ps 



d 2 r [(V??) 2 + (V6») 5 



(4) 



with rigidity p s = h 2 n/m*. The vortex elastic energies 
and vortex interactions in logarithmic approximation can 
be calculated from ^ similarly to the usual case 

The energy of a single vortex is E^f = np s (k 2 + 
m 2 )\n(R/a), where a = h/(2m* p) 1 / 2 is the core radius 
(see below) and R 3> a is the radius of the polariton exci- 
tation spot. Clearly, the elastic energy of the half- vortex 
is twice smaller than the energy of a usual vortex [e.g., of 
the (0, 1) vortex]. The energy E^' of the pair of vortices 
(fci,mi) and (&2, 777.2) separated by distance r, such that 
a <C r <C R, is 

= 7rp s [(fci + fc 2 ) 2 + (mi + 777 2 ) 2 ] \n(R/a) 

+ 2irp s (k 1 k 2 + mim 2 ) ln(a/r). (5) 



and left-circular for sgn(fcm) = — 1. In particular, there 
is attraction between the right half-vortex and the right 
anti-half-vortex, (1/2,1/2) and (-1/2,-1/2), as well as 
between the left ones, (1/2,-1/2) and (-1/2,1/2). On 
the same time, the right half- vortices do not interact with 
the left half-vortices, and, for example, the energy of the 
pair (1/2,1/2) and (1/2,-1/2) simply equals to the en- 
ergy of a single (1,0) vortex. [13j 

At a finite temperature equal numbers of left and right 
half- vortex pairs are excited. Since the subsystems of left 
and right vortices do not interact, they evolve indepen- 
dently with increasing temperature and are subject to 
the Kosterlitz-Thouless transition at T c — (jr/4)p s - This 
temperature is twice smaller than the critical tempera- 
ture for spin-less bosons due to the double reduction of 
the single half-vortex energy. Clearly, this estimation of 
T c does not take into account the depletion of the con- 
densate, but it shows the necessity to allow for the po- 
larization degree of freedom for any realistic calculations 
of the transition temperature. The critical temperature 
is expected to be modified also by non-parabolicity of 
polariton kinetic energy, as well as by the account for 
the longitudinal-transverse splitting of polariton branch 
that lead to the coupling between left and right vortices. 
Note also that the symmetry between left and right vor- 
tex subsystems is broken in applied magnetic field. 

Polarization texture of half-vortex core. — The chemi- 
cal potential p is found experimentally as a blue-shift 
of polariton luminescence line due to the polariton con- 
densation, and, typically, p < 1 meV [3, hj. Due to a 
very small value of the polariton effective mass to* rang- 
ing from 10 -4 to 10~ 5 of the free electron mass [151 ] . the 
size of half-vortex core is quite big, a > 1 pm. There- 
fore, the polarization texture within the core region of 
a half- vortex can be observed by means of near-field lu- 
minescence spectroscopy. This texture is studied in this 
subsection. 

In what follows, only the basic half-vortices with |fc| = 
|m| = 1/2 will be considered. In this case, the order 
parameter can be written in cylindrical coordinates as 

tfhv = [A(0/(r/a) - iB(^) 5 (r/o)] , (6) 

where the azimuthal dependencies are given by 



A(0) = e im ^{cos(fc0),sin(fc</>)}, 



(7a) 



B(4>) = sgn(fcm)e im *{sin(fc(/.), - cos(fc0)}, (7b) 

and the radial functions f(r/a) and g(r/a) are found 
from minimization of the vortex energy, i.e., from the 
equation 8H/8%p = 0. Using the relations 



It is seen that only half-vortices with the same sign of 
the product km interact. This sign, as we see below, 
defines the sign of circular polarization in the center of 
the vortex core — it is right-circular for sgn(fc777,) = +1 



A»(0)=iB"(0) = --(A + iB) ! 



A — iB — {1, sgn(fcm)i}, 



(8a) 



(8b) 
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FIG. 2: Showing the half- vortex radial functions f(r/a) 
(dashed) and g(r/a) (solid) for three values of the interac- 
tion parameter 7. 



one obtains (£ = r/a) 

f" + \f 2|2 (/ - 9) + [l-f- 19 2 }f = 0, 



(9a) 



1 

2f 2 



o(9 -f) + [1-9 -7/19 = 0, (9b) 



where 7 = (<7 + C/i)/(t/ - E7i). 

Equations (J5^,b) are symmetric with respect to inter- 
change of radial functions / and g. The half- vortex, how- 
ever, is described by an asymmetric solution. Since the 
polarization should become linear at large distances, one 
has /(oo) = 1 and 3(00) = 0. This way at »• > a the 
solution ([51) transforms into ^ with ry = kip and 9 = m<f>. 
On the other hand, in the half-vortex center one has to 
demand /(0) = g(0) in order to remove the divergences 
produced by terms (/ — g)/2£ 2 in Eqs. ([9]). Taking into 
account the relation (JSb) it is clear that the polariza- 
tion is fully circular in the core center (r = 0), and the 
sign of circular polarization is given by sgn(fcm), as was 
mentioned above. 

Equations ([9]) have simple solutions in the particular 
case U\ — Uq/2, i.e., for 7 = 3. Namely, /(£) = [1 + 
h(Cl\/2 and g(£) = [1 - h{£)}/2, wh ere h(£) is the radial 
function of a usual vortex in spin- less condensate 
Since h(£) is monotonously decaying from 1 at £ — > 00 to 
at £ = 0, the polariton concentration n(r) — f 2 (r/a) + 
g 2 (r/2) is decreasing from its unperturbed value n = 
n(oo) to n/2 in the half- vortex center 



In general, for realistic values of 7 [17[ the behavior of 
/ and g remains close to the above case. The numeri- 
cal solutions of Eqs. ([9]) are shown in Fig. [2j Note that 
while the asymptotic behavior of /(£) is 7-independent 
for large £, /(£) ~ 1— (4£ 2 ) -1 , the second radial function 
behaves asymptotically as g(£) ~ [2(7 — 1)£ 2 ] -1 . There- 
fore, the effective size of the half- vortex core grows with 
decreasing 7 and diverges at 7 — > 1, i.e., when Ux — > 0. 



This reflects the fact that polaritons are not expected 
to exhibit the superfluid transition for U\ = 0. In this 
case the polariton system has 0(4) symmetry and the or- 
der is destroyed at any finite temperature according the 
nonlinear cr-model analysis. 

Concerning the polarization texture of the half- vortices 
given by Eqs. ([6]) and (0), it should be noted that dif- 
ferent half-vortices can be transformed to each other by 
applying two symmetry operations, (i) the 2D inversion 
(x — > x, y — > —y), and (ii) the time inversion or complex 
conjugation. Each of these operations used separately 
transforms the right (1/2,1/2) half- vortex into one of 
two left half-vortices, while the subsequent application of 
both operations yields the (—1/2, —1/2) anti-half-vortex. 
Clearly, these symmetry operations leave the radial func- 
tions f(r/a) and g(r/a) unchanged. 

The half-vortices possess two anisotropic polarization 
textures shown in Fig. [3ja,b). These figures rely on the 
usual representation of the time-dependent electric field 
as E(£) oc Rc{xpe~ luJt } , where Li is given by the bare 
frequency luq of the lower-branch polariton in microcav- 
ity blue-shifted by the chemical potential fi due to the 
polariton-polariton repulsion, u> — ujq + /i. It is seen that 
when one approaches the core center linear polarizations 
convert into elliptical ones. The circular polarization de- 
gree is given by 



sgn(fcm) 



2/ff 
f 2 +9 2 



(10) 



and |pcirc| increases from to 1 with decreasing r. In 
spite of the dependence of the directions of main axes 
of each elliptical polarization on the azimuthal angle 0, 
all elliptical polarizations meet in phase at r = and 
transform into the same circular polarization. Note also 
that each image in Fig.[3]represents two half- vortices, one 
with the clockwise rotation of polarization vector and the 
other with the counterclockwise rotation. 

The above analysis of half- vortex properties is based 
on minimization of the polariton energy subject to spe- 
cific boundary conditions, which is valid for the case of 
polariton condensate in quasi-thermal equilibrium. Such 
thermalized polariton condensates have being recently 
observed 0, [lj| • In these conditions the half- vortex pairs 
can both appear spontaneously and be artificially ex- 
cited. In particular, the half-vortex and anti-half-vortex 
pair can be obtained by shining the uniform linearly po- 
larized condensate with two circularly polarized pulses 
having appropriate intensities, spot radii, and separa- 
tion. Another option consists in the excitation of a (0, 1) 
vortex using the optical parametric oscillator setup [191 ] . 
This vortex involves only phase change and is topologi- 
cally equivalent to the pairs of (1/2, 1/2) and (—1/2, 1/2) 
half- vortices, or to the pair of (-1/2, 1/2) and (1/2, 1/2) 
ones, so it can decay in one of these pairs. 

Finally, it should be noted that similar half-vortices 
with two topological charges can be present in mixtures of 



4 




FIG. 3: Showing two distinct core textures of basic half- 
vortices. The case (a) is realized for the right (1/2, 1/2) and 
the left (1/2, —1/2) half-vortex, while the case (b) — for right 
(-1/2,-1/2) and the left (-1/2,1/2) half-vortex. Arrows 
indicate the instant polarization at two values of r = const 
(dashed circles) and in the core center r = 0. The polarization 
vectors change in time following red solid lines. 

two hyperfine spin states of 87 Rb atomic condensates [2(| , 
provided the order parameter manifold is given by two 
angles as in Fig. [T] In this case the half-vortices should 
also control the Kosterlitz-Thouless transition observed 
in these condensates (2lj . 

Conclusions. — Polarization vortices in linearly polar- 
ized polariton condensate have been analyzed. It has 
been shown that four half- vortices with the direction and 
phase winding numbers equal to ±1/2 have the smallest 
energy. These half- vortices form two subsystems exhibit- 
ing independent Kosterlitz-Thouless transitions at the 
same temperature. The condensate polarization becomes 
elliptical within the core of a half- vortex and converts to 
fully circular in the core center. 
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